Signals of Inert Doublet Dark Matter in Neutrino Telescopes 
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One of the simplest extensions of the Standard Model that explains the observed abundance of 
dark matter is the inert doublet model. In this theory a discrete symmetry ensures that the neutral 
component of an additional electroweak doublet scalar is stable, and constitutes a dark matter 
candidate. As massive bodies such as the Sun and Earth move through the dark matter halo, dark 
matter particles can become gravitationally trapped in their cores. Annihilations of these particles 
result in neutrinos, which can potentially be observed with neutrino telescopes. We calculate the 
neutrino detection rate at these experiments from inert doublet dark matter annihilations in the 
cores of the Sun and the Earth. 



I. INTRODUCTION 

One of the most pressing questions currently in physics 
is the nature of dark matter. It is known that most of 
the matter in the universe is comprised of non-luminous 
(dark), neutral matter. Moreover, it is also known 
that none of the standard model particles can be the 
dark matter. Hence, to explain dark matter, we must 
introduce new physics beyond the standard model. The 
most appealing scenario from a cosmological standpoint 
is that of a weakly interacting massive particle (WIMP). 
The WIMP is theoretically attractive because any weakly 
interacting particle with a mass in the TeV scale will give 
approximately the correct relic abundance of dark mat- 
ter. Many extensions of the standard model contain such 
dark matter candidates, for example, supersymmetry [Q, 
extra dimensions [| and little Higgs Q . 

Dark matter is distributed across the galaxy in the 
form of a halo. As the dark matter particles interact with 
nuclei in large objects such as the Sun or Earth, a certain 
fraction of the scattered dark matter particles will not 
have enough kinetic energy to escape the gravitational 
pull of the object. Thus, eventually a build up of dark 
matter will occur in the cores of massive objects. The 
dark matter particles will then annihilate, producing 
standard model particles, which eventually decay to 
neutrinos. If the neutrinos originate from the cascade of 
a TeV scale object such as a WIMP, they will typically 
have an energy of a few hundred GeV. This energy range 
is detectable in neutrino telescopes such as IceCube. 
This idea has been investigated before in the context of 
supersymmetric dark matter [^| and Kaluza-Klein dark 
matter candidates 

In this paper, we consider the possible signal in 
neutrino telescopes from dark matter annihilations in the 
Sun and Earth for the inert doublet model (IDM) . In this 
model, the Higgs sector of the standard model (SM) is 
extended to include an additional SU(2) doublet that 
does not acquire a vacuum expectation value (vcv) . If an 
additional Zi symmetry is imposed, the neutral particle 
in the additional doublet becomes stable and thus, a 
dark matter candidate. The possible neutrino signal from 
the inert doublet model is especially interesting because 



this type of dark matter arises in extensions of the SM 
motivated by the little hierarchy problem, such as the 
Left-Right Twin Higgs model 0. 

The format of the paper is as follows: In Section II we 
discuss the basics of the IDM and list some constraints. 
In section III we discuss the dynamics of capture of dark 
matter particles by massive objects and calculate the 
capture and annihilation rates for the Sun and Earth 
for the inert doublet model. In Section IV, we discuss 
the determination of the detector rate. In section V we 
present our results, and in section VI we conclude. 



II. THE INERT DOUBLET MODEL 

The inert doublet model consists of the standard model 
with an additional electroweak doublet scalar, H2, with 
the same quantum numbers as the SM Higgs. This 
doublet is odd under a Z2 symmetry, while all other fields 
are even, and does not acquire a vev. Hence, H2 couples 
only to SM gauge bosons, the SM Higgs and itself. Since 
H% is odd under the Z2 symmetry, this ensures that the 
lightest component of H2 is stable. Thus, if the neutral 
scalar is lighter than the charged scalar in H2, it is a good 
dark matter candidate |7], || . 

The general pattern of symmetry breaking in a two 
Higgs doublet model was first investigated by Deshpande 
and Ma ||. Recently, Barbieri et al. have shown that 
in the case of an unbroken Z2 and appropriately chosen 
splittings the model passes electroweak precision tests 
with a heavy Higgs, thus solving the little hierarchy 
problem . They refer to this case as the inert doublet 
model since the new Higgs doublet cannot couple to 
fermions. The IDM has been shown to give the correct 
relic abundance for dark matter jyj. The direct detec- 
tion rates of dark matter in this model were examined 
by |p^| . It has also been analysed for indirect signals of 
dark matter via photons 1 13 and monochromatic photon 
production at the Galactic center jjS. The LEP II 
analysis for neutralino pair production has recently been 
translated into constraints on the IDM |15|. 



A. Parametrization of the Potential 

The most general potential consistent with the Z2 
symmetry is 

V = filH!] 3 + i4\H 2 \* + 
+ A 2 |i/ 2 | 4 + A 3 | J ff 1 | 2 | J ff 2 | 2 

+ X i \HlH 2 \ 2 + ^-{(HlH 2 f + h.c]. (1) 
Expanding the potential above in unitary gauge, 

Hl = {(v + lyy/2) (2) 

we obtain, 



with WMAP at the 3er level: a low mass region, where 
ms < 100 GeV, and a high mass region, where 500 GeV 
< m Ho < 2 TeV. 



2. From direct detection bounds 

The direct detection bounds put a constraint on the 
value of A5. As A5 — > 0, the neutral scalar S and 
the pseudo-scalar A become degenerate. In this case, 
these particles can scatter off matter via a Z exchange, 
giving cross-sections eight to nine order of magnitudes 
larger than present direct detection bounds This 
constraint can be avoided if the mass splitting between 
S and A is higher than the kinetic energy of dark-matter 
in the halo, or 82 ~ a few hundred MeV. 



3. From V z 



rrr 



= 2w 2 Ai = -1u\ 

= A + \(h + A 4 + X 5 )v 2 

= M2 + ^( A 3 + A 4 - A 5 )u 2 

2 1 \ 2 
= M 2 + 2 X 3 V ■ 



(4) 



This leaves us with seven independent real parameters: 
Hi, Ai, A2, A3 , A4, and A5. Fixing the Z mass fixes v, 
while fixing the SM Higgs mass fixes Hi and Ai . Following 
JiT| we define A^ = (A3 + A4 + A5). Fixing the scalar mass 
and the mass splittings between the scalar, pseudoscalar, 
and charged particles as ms, 81 = m H ± — ms and 82 — 
mA — ms, fixes /i2, A3, and A5. This leaves us with the 
new parameter set: mz, mu, A2, Xl, 81, 82, and ms- 
The Z mass is fixed by LEP, while we take A2 = 0.1 
p3| . This parameter does not affect the neutrino flux 
rate calculation directly. 



B. Constraints 

There are a number of constraints that serve to limit 
our parameter space. We discuss each of them in detail 
below. 



1. From Qh 2 

Recent measurements from WMAP have led to a 
precise determination of the amount of dark matter in 
the Universe: VL DM h 2 = 0.112 ±0.009 jflfJ. The program 
micrOMEGAs |7| uses a calcHEP Jl§ model file to 
solve the Boltzman equation numerically to find the 
relic density. We used this to exclude a large portion 
of the parameter space. Two regions are consistent 



Data from LEP has put tight constraints on the width 
of the Z boson, such that new decay channels are very 
restricted. This puts lower bounds on the scalar mass 
and the mass splittings. 



ms + mA > mz => ms > 



2m H ± > mz =>■ ms > 



m z - 8 2 



mz — 2Si 



4- From model stability 

In order to ensure stability of the model, we require 
Al,A3 > — 2 v / AiA2- Since we consider the scalar, S, to 
be the dark matter candidate, Xl < A3, and the first 
condition is sufficient. Substituting the SM Higgs mass 

for Ai yields the condition: X L > -^y^f-. For a 
modest choice of A2 and m^, this excludes almost all 
negative Al couplings. 



III. DARK MATTER CAPTURE AND 
ANNIHILATION 

Dark matter particles (DMPs) accumulate in the 
massive bodies from the Galactic halo, and are depleted 
by annihilations. Given a long enough time, this process 
can come into equilibrium. The differential equation 
governing the number of DMPs in the Sun (or Earth) 
is 



N=C- C A N 2 (5) 

Here C is the capture rate from the halo, and Ca = 
(<jav)V2/V 2 . (<tav) is the total cross-section times the 
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relative velocity in the limit v — > 0. Vj are the effective 
volumes for the Sun and Earth, given by 



Vj = (3m 2 pl T/(2jm s p)) 3 / 2 



(6) 



and rripi is the Planck mass, T is the temperature and 
p the core density of the massive body. Equation (|^) can 
be solved to obtain the annihilation rate of DMPs 



T A = -Ctanli (i/r) 



(7) 



If the time scale for equilibriation (r = \[UCa) is much 
smaller than the age of the solar system (t) , then Ta is 
simply Since this is the scenario of a maximal signal, 
this indicates that the capture rate is the dominant factor 
in determining the magnitude of the signal. 



A. Capture Rate 

The capture rate depends on the DMP scattering 
cross-section from the nuclei. The calculation of cross- 
section proceeds in three steps. First we calculate the 
partonic cross-section with quarks and gluons. We next 
translate this into the interaction with nucleons, by 
taking the quark/gluon matrix elements in the nucleonic 
state. Finally, we evaluate the nucleon operator matrix 
elements in the nucleus [p0| . Since these particles are 
non-relativistic, the cross-section calculation is greatly 
simplified. The DMP undergoes scalar interactions which 
coherently add in the nucleus, and hence it couples to 
the mass of the nucleus. The elastic cross-section at zero 
momentum transfer is [11(1 



? ? 

m S m N 



47r(ms + rriNi) 2 V^sm^ 



A, 



f 2 m 



(8) 



Here, mjv is the mass of the nucleus and / ~ 0.3 is the 
nucleonic matrix element 21 , defined by 



(N\J2m q qq\N) = fm N {N\N). 



(9) 



At a finite momentum transfer, the particle doesn't 
see the entire nucleus, and hence the cross-section is 
suppressed by a form-factor (Fi(ms)). Other factors 
that influence the capture rate are the elemental abun- 
dance fi, distribution (pi and the kinematic suppression, 
^(ms/mjvj p2]. The capture rate is given by 



C = 



Po.3 



x fi4> 



(m s /GeV)v 2 7Q 
S(m s /m Ni ) 



10- 4 pb 



(m N JGeV) 



(10) 



where c = 4.8 x 10 24 s" 1 for the Sun and c = 4.8 x 10 15 
s _1 for Earth, po.3 is the local halo mass density in units 
of 0.3 GeV cm -3 , and U270 is the dark-matter velocity 
dispersion in terms of 270 km s _1 . The quantities fi and 
(j>i for the Sun and Earth are given in p3fl The sum is 
over all elements in the Sun or Earth. The kinematic 
suppression factor can be parametrized as following 



S(x) 
A 



[Ayii + Ai)}* 



2(x- 1) 



(11) 
(12) 



The mean escape velocity, (v esc ) is obtained by numer- 
ical fitting to the exact kinematic suppression. For the 
Sun, (v esc ) = 1156 km s _1 , and for the Earth (v esc ) — 
13.2 km s _1 . This factor arises from the fact that in 
order to be captured, the scalar has to scatter off a 
nucleus with a velocity less than the escape velocity 
in that object. If the mass of the nucleus is equal 
to the mass of the scalar, then it can potentially lose 
all its kinetic energy in an elastic collision. It can be 
seen from the parametrization above that the capture is 
suppressed more if the mass of the dark matter particles 
is very different from the mass of the nucleus. The 
form factor suppression can be modeled numerically for 
simplicity. For the Earth, the only appreciable effect is 
the suppression from iron, so Fj(ms) ~ 1 for all other 
elements. For iron, 



Fpe 



1 - 0.26 



.4 



A+l 



(13) 



For the Sun, the fitting is slightly more complicated. 
The suppression is given by 



F i (m s ) = ^r / + (l-^r / )exp 



inf\ 



lpg(TOs) 

log(m 4 c ) 



(14) 



The quantities i 7 !/"^, mt. and on are numerical fitting 
parameters taken from p3|. 



B. Annihilation 

These accumulated DMPs can annihilate into Stan- 
dard Model particles. Subsequently, these particles 
can decay into energetic muon neutrinos, and hence be 
detected in astrophysical neutrino detectors. The scalar 
S usually annihilates into a two-body state, and hence 
the typical energy of the neutrino is ~ jins- 

The overall annihilation rate in the Sun or Earth 
is essentially set by the capture rate (assuming near 
maximal signal). We calculate the dark matter annihi- 
lation branching fractions into various Standard Model 
particles. We can then use known parametrizations for 
neutrino spectra from decay of these particles in the Sun 
or Earth. The branching fractions for a specific set of 
parameters is shown in Figure |l|. 
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FIG. 2: Second moment of neutrino spectra from the Earth (left) and the Sun (center). Also shown is the second moment of 
anti-neutrino spectrum from the Sun (right). The Earth anti-neutrino spectrum is identical to the neutrino spectrum. 



IV. NEUTRINO SPECTRA AND DETECTOR 
RATE 

The neutrino spectra from annihilation events de- 
pend upon solar parameters, neutrino physics, quark 
hadronization etc. Once we know the neutrino 
spectra from each of the decay products, the overall 
spectrum is given by convolving these spectra with the 
branching fractions. For the purposes of measurement 
in the neutrino telescopes, the relevant quantity is the 
detector rate. The technique for inferring the existence of 
the neutrino is observation of a muon, which is produced 
by a charged-current interaction. The cross-section for 
this process is proportional to the energy of the neutrino, 
and the range of the subsequent muon is also proportional 
to its energy. Thus, we are interested in the second 
moment of the neutrino spectrum for the detector rate 
as a function of the injection energy (E in ) p5|, |26|. 



(N Z 2 ) F ,i{E in ) 



dN 
~dE 



FA 



E 2 

(E u ,E in )-g-dE v (15) 
E im 



The spectra from the Sun are more complicated than 
the Earth spectra. For the Earth, we consider the 
neutrino spectrum in the rest frame of the decaying 
particle, and then boost it for a particle with an energy 



Ei n . If the injected particle is a b- or c-quark, we also 
have to take hadronization into account. The quark 
loses energy as it hadronizes, so the injected energy is 
a fraction of Ei n . In the Solar case, we also have to 
consider the stopping of heavy hadrons. The core of the 
Sun is dense enough to slow b- and c-quarks further after 
hadronization. Another effect in the Sun is neutrino 
stopping and absorption. Neutrinos lose energy via 
neutral-current interactions and can be absorbed through 
charged-current interactions in the Sun. Stopping and 
absorption co-efficients turn out to be different for neu- 
trinos and anti-neutrinos. Hence, we expect a distinct 
spectrum for neutrinos and anti-neutrinos from the Sun. 
From the Earth, these spectra are identical. We produce 
the (Nz 2 ) plots in Figure |^. Detailed numerical fits 
for these functions can be found in p3| . The combined 
detector rate is given by 



detect 



£)( 



GcV 



^2a l b l ^B F {Nz 2 ) F ^{m s ) 



(16) 



The sum over i is over neutrino and anti-neutrino 
states. The are the neutrino scattering co-efficients, 
a v = 6.8, a p = 3.1. The bi are the muon range 
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FIG. 3: Detector rate from annihilations in the Earth (left) and the Sun (right) in the low mass region for various values of 
\l- Different plots are for different Higgs masses and mass splittings (81 and 82)- The solid sections in the lines represent the 
mass range consistent with the correct relic abundance. 
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FIG. 4: Detector rate from annihilations in the Earth (left) and the Sun (right) in the high mass region for various values of Az,. 
Different plots are for different mass splittings (Si and 82). The solid sections in the lines represent the mass range consistent 
with the correct relic abundance. 
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co-efficients, b v = 0.51 and 6 P = 0.67. Bp is the 
annihilation branching fraction of the DMP for channel 
F. 

For the sun, the constant c = 2.54 x 10~ 23 km~ 2 yr" 1 . 
The expression for the detector rate for the earth is scaled 
by the square of ratio of the earth-sun distance to the 
earth radius. Thus, c = 1.42 x 10~ 14 km~ 2 yr _1 for the 
earth. 

The background for this process arises mainly from 
atmospheric neutrinos J27J]. For the case of the Sun, 
we can reduce the background by including events only 
within a narrow angular cone along the line of sight 
from the detector to the Sun. There is also some 
background from the solar atmospheric neutrinos p8| , 
where neutrinos are produced in the solar atmosphere by 
cosmic rays. 



V. RESULTS 

We have calculated the expected rate of detection in 
neutrino telescopes from these annihilations in the IDM. 
We restrict the parameter space to that allowed by relic 
abundance constraints by WMAP and other constraints 
discussed above. The IDM allows the Higgs mass to 
be higher than the standard electro-weak precision test 
constraints jlQfl , so we also analyze the signal for a 200 
GeV SM Higgs. We do the calculation for different sets 
of mass splittings, S\ and 62- The plots are produced as a 
function of ms for various values of Al- Wc find that the 
results are qualitatively different for the low mass (ms < 
100 GeV) and the high mass (500 GeV < m s < 2 TeV) 
regions, hence we present them separately. 

We plot the number of events observed per year in a 
detector with an effective area of 1 km 2 , such as IceCube. 
The solid sections on the lines indicate the mass range 
which gives the correct relic abundance. 



A. Low Mass Region 

In the low mass region, the signals from the Earth and 
the Sun are comparable. The maximal signal for Earth is 
obtained around m$ — 50 GeV. The various peaks in the 
Earth plots are explained by kinematics. When the DM 
mass is close to a nucleus mass there is no kinematic 
suppression (Eq. jll|)). As the DM mass increases, 
the cross section for DM capture decreases, leading to 
a reduction in the rate. In the Solar case, there is no 
kinematic resonance, therefore the Earth signal exceeds 
the Solar signal at certain to 5 values. 

The branching fractions in the low mass region (below 



the W threshold) , where annihilation products are quarks 
(b and c) and leptons (t), are independent of A^. The 
detector rate depends on Al only through (&av) and the 
capture rate (C) in this region. 

We see ~ 100 events from the Earth with a 200 GeV 
Higgs and tens of events with a 120 GeV Higgs (Figure 
|3). We expect to see a few hundred events per year from 
the Sun in this parameter range. 



B. High Mass Region 

In the high mass region, the signal from Earth is 
suppressed by several orders of magnitude as compared 
to the signal from the Sun. This is because the annihi- 
lation rate is not maximal (the capture and annihilation 
processes have not come into equilibrium) for the Earth, 
but it is for the Sun. The detector rate decreases with 
the mass of the dark matter particle, for both Earth and 
the Sun, due to the cross-sectional dependence on ms 
and kinematic suppression. The branching fractions for 
the high mass particles have a non-trivial dependence 
on Xl. Consequently, the qualitative features of the 
detector rate are influenced by the branching fractions 
and annihilation rate which, in turn, depend on the 
specific value of Al. 

In this mass range, the signal from the Earth is too 
low to be observed. The signal from the Sun is more 
promising, and we expect to see a few events per year. 



VI. CONCLUSION 

Neutrino telescopes provide a very interesting mech- 
anism for dark-matter detection. With the increased 
sensitivity of IceCube, these signals can be used to test 
large ranges of parameter space in various models. We 
find a promising signal in the inert doublet model in two 
distinct parameter ranges. We get a few hundred events 
per year in the low mass region, and a few events from 
the Sun in the high mass region. For detection, a careful 
analysis of signal over background may be required. 
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